Preprint typeset in JHEP style - HYPER VERSION 



arXiv:1009.0218 



Constraining warm dark matter with cosmic shear 
power spectra 



Katarina Markovic a fe , Sarah Bridle c , Anze Slosar d , Jochen Weller a b e 

a University Observatory Munich, Ludwig- Maximilian University, Scheinerstr. 1, 81679 
Munich, Germany 

b Excellence Cluster Universe, Boltzmannstr. 2, 85748 Garching, Germany 

c Department of Physics and Astronomy, University College London, Cower Street, 

London WC1E 6BT, United Kingdom 

d Building 510 A, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 
e Max- Planck- Institut for Extraterrestrial Physics, Giessenbachstr. , 85748 Garching, 
Germany 

E-mail: markovic(at)usm. lmu.de, sarah.bridle(at)ucl. ac.uk, 
anze (at ) bnl . gov , j ochen . weller (at ) usm . lmu . de 

Abstract: We investigate potential constraints from cosmic shear on the dark matter 
particle mass, assuming all dark matter is made up of light thermal relic particles. Given 
the theoretical uncertainties involved in making cosmological predictions in such warm 
dark matter scenarios we use analytical fits to linear warm dark matter power spectra 
and compare (i) the halo model using a mass function evaluated from these linear power 
spectra and (ii) an analytical fit to the non-linear evolution of the linear power spectra. 
We optimistically ignore the competing effect of baryons for this work. We find approach 
(ii) to be conservative compared to approach (i). We evaluate cosmological constraints 
using these methods, marginalising over four other cosmological parameters. Using the 
more conservative method we find that a Euclid-like weak lensing survey together with 
constraints from the Planck cosmic microwave background mission primary anisotropies 
could achieve a lower limit on the particle mass of 2.5 keV. 
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1. Introduction 



In the second half of the 20th century, two competing theories for the growth of cosmolog- 
ical structure were proposed. In the cold dark matter (CDM) paradigm (Peebles [1982], 
Blumenthal et al. [1984], Peebles [1984], Davis et al. [1985]), dark matter possesses a neg- 
ligible initial velocity dispersion and structure formation proceeds in a strictly hierarchical 
bottom-up manner: small haloes of dark matter form and later merge into larger ones, 
which in turn merge into even larger haloes. In these virialised dark matter structures 
the baryons condense and form luminous objects in the Universe. In the hot dark matter 
(HDM) paradigm (Zel'Dovich [1970], Bond et al. [1980], Bond and Szalay [1983], Centrella 
et al. [1988]) on the other hand, the dark matter is ultra- relativistic in the very early stages 
of the Universe, erasing all structure on small scales. In these models, the most massive 
structures form first, producing "Zeldovich pancakes", that later produce smaller objects 
by fragmentation in a top-down manner. An example of such an extremely energetic dark 
matter particle is a massive active neutrino. 

By the end of the twentieth century it was clear that the hot dark matter paradigm 
cannot describe the measurements of the cosmic microwave background and the clustering 
of galaxies and that structure formation in the Universe is, at least overall, hierarchical 
(Komatsu et al. [2010], Cole et al. [2005], Tegmark et al. [2004], Seljak et al. [2006a]). 

However, there exist a number of observations that may indicate a lack of small scale 
perturbations in the initial density field, relative to that predicted in the standard ACDM 
paradigm. For example, it has long been known that CDM theory predicts many more 
small mass haloes than the number of dwarf galaxies that we see around the Milky Way 
(Diemand et al. [2007] and references therein). Similarly, cuspy galactic cores indicated 
in some observations are inconsistent with predictions of the CDM (Moore [1994], Simon 
et al. [2005]). Moreover, the angular momenta of dark matter haloes are considerably lower 
than those observed in spiral galaxies (Sommer-Larsen and Dolgov [2001], Chen and Jing 
[2002], Zavala et al. [2008]). There is also some discrepancy between the distribution of 
sizes of mini- voids in the local Universe and CDM predictions (Tikhonov et al. [2009]). 

These discrepancies might be resolved by accounting for certain astrophysical pro- 
cesses. Supernova feedback can extinguish star formation and further baryonic effects can 
also affect the properties of the dark matter density distribution in centres of haloes. How- 
ever, a suppression of the primordial matter power spectrum on small scales is an attractive 
alternative. This is most easily achieved by giving dark matter some small initial velocity 
dispersion: not enough to break the very successful hierarchical structure formation, but 
enough to make a difference on small scales. Such models go under the name of warm dark 
matter (WDM) (Bode et al. [2001], Avila-Reese et al. [2001]). 

In warm dark matter models, dark matter particles free-streamed for a short period 
in the early Universe, before becoming non-relativistic. During this period, fluctuations 
on small scales were heavily suppressed and hence WDM models are characterised with a 
free-streaming wavenumber kf s , beyond which the fluctuations in the power spectrum of 
initial fluctuations are exponentially suppressed. This suppression is the main observational 
smoking gun of WDM models. For the models discussed in this paper, Af s = 2ir /k{ s is of 



- 1 - 



the order of a Megaparsec. 

Several microscopic models for warm dark matter have been proposed. The most 
common models contain sterile neutrinos (Dodelson and Widrow [1994], Fuller et al. [2003], 
Asaka et al. [2005], Abazajian [2006], Boyarsky et al. [2006b], Petraki and Kusenko [2008], 
Laine and Shaposhnikov [2008], Kusenko [2009], Hamann et al. [2010]) or gravitinos (Bond 
et al. [1982], Borgani et al. [1996], Fujii and Yanagida [2002], Cembranos et al. [2005], 
Steffen [2006], Takahashi [2008]) as dark matter particles. 

In sterile neutrino models, the standard model of particle physics is extended with the 
addition of one or more singlet right-handed sterile neutrinos. The standard model active 
neutrinos then often receive the observed small masses through a see-saw mechanism. This 
implies that sterile neutrinos must be rather heavy, but the lightest of them naturally has 
a mass in the keV region, which makes it a suitable warm dark matter candidate. The 
simplest model of sterile neutrinos as WDM candidates assumes that these particles were 
produced at the same time as active neutrinos, but they never thermalised and were thus 
produced with a much reduced abundance due to their weak coupling. In the second class of 
models the gravitino appears as the supersymmetric partner of the graviton in supergravity 
models. If it has a mass in the keV range, it is a suitable WDM candidate. It belongs 
to a more general class of thermalised warm dark matter candidates. It is assumed that 
this class of particles achieved a full thermal equilibrium, but at an earlier stage, when the 
number of degrees of freedom was much higher and hence their relative temperature with 
respect to the CMB is much reduced. 

When describing WDM structure formation there are three relevant variables that 
describe a particular model of WDM: the temperature Twdm of WDM particles, their 
mass w-wdm an d their energy density wwdm- However, the equations of the evolution 
of perturbations can be re-parametrised as a two-parameter family of models, where one 
typical specifies wwdm = S^wdm^ 2 (with h = Hq/10Q km s _1 Mpc _1 denoting the reduced 
Hubble's constant and fiwDM for the fractional energy density of WDM in units of critical 
density) and the ratio ttt-wdm / ?wdm • This allows one to make the following one-to-one 
relation between the mass of the thermalised WDM particle uiwdm, such as the gravitino, 
and the mass of the simplest sterile neutrino m us , such that the two models have an identical 
impact on cosmology (Viel et al. [2005]) 



The difference comes from the fact that in the gravitino case the particle is fully thermalised, 
the number of effective degrees of freedom being determined by mass and energy density of 
dark matter, while in the simplest sterile neutrino case the number of degrees of freedom 
is fixed, while abundance is determined by mass and energy density of dark matter. In the 
rest of this paper we work with fully thermalised WDM particles. 

The above WDM models have been constrained using several observational techniques 
(e.g. Polisensky and Ricotti [2010], Boyanovsky et al. [2008a] etc.). The most popular is 
the observation of the flux power spectrum of the Lyman-a forest, which can indirectly 
measure the fluctuations in the dark matter density on scales between ~ 100 kpc and 
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~ 10 Mpc. The limits are ttt-wdm > 4 keV or equivalently m us > 28 keV at 95% confidence 
limit (Viel et al. [2008], see also Viel et al. [2005], Seljak et al. [2006b]). For the simplest 
sterile neutrino model, these lower limits are at odds with the upper limits derived from 
X-ray observations, which come from the lack of observed diffuse X-ray background from 
sterile neutrino annihilation and set the limit m vs < 1.8 keV at 95% confidence limits 
(Boyarsky et al. [2006a]). We stress that these results rule out the simplest sterile neutrino 
models. However, there exist theoretical means of evading small-scale power constraints 
(see e.g. Boyarsky et al. [2009] and references therein). 

The distribution of matter can be observed using the bending of light by gravity, 
termed gravitational lensing. This causes a distortion in the apparent shapes of distant 
galaxies, which is described to first order by a shear. The effect is generally small but can be 
measured by statistically averaging many observations together, termed cosmic shear. The 
two point correlation function of the shears of pairs of galaxies is a probe of the underlying 
matter power spectrum and thus of the nature of dark matter. 

The latest cosmic shear measurements (Fu et al. [2008]; Massey [2009], Schrabback 
et al. [2009]) have sufficient precision to just constrain one parameter for the amplitude 
and one for the shape of the matter power spectrum, usually written as the linear theory 
rms fluctuation in 8 /i" 1 Mpc spheres, a% and the matter density of the Universe, f2 m . 
Cosmic shear has been identified as the method with the most potential to solve one of the 
major problems in cosmology: the nature of dark energy (Albrecht et al. [2005]; Peacock 
et al. [2006]). There are therefore several upcoming large cosmic shear surveys both from 
the ground (Pan-STARRS 1 , the Dark Energy Survey 2 , HyperSuprimeCam 3 and the Large 
Synoptic Survey Telescope 4 ) and space (Euclid 5 and WFIRST 6 ), with additional surveys 
proposed. 

Non-gravitational interaction between dark matter and baryonic matter is extremely 
weak or non-existent. Therefore, gravitational lensing is one of the most powerful probes 
of the nature of dark matter. Because it can detect the presence of dark matter directly, 
it provides some strong evidence for its existence (Clowe et al. [2004]). Measuring grav- 
itational lensing also enables us to map-out the 2 and even 3-dimensional distribution of 
dark matter (Kaiser and Squires [1993], Hu [1999] and Taylor et al. [2004]). We can also 
calculate the statistical properties of the matter distribution, which allow us to infer the 
matter density and the root-mean-squared amplitude of matter density fluctuations (Hoek- 
stra et al. [2002]). Since dark matter particle properties may directly influence the matter 
distribution, these measurements can additionally constrain the model for dark matter 
particle origin and nature (e.g. Schaefer et al. [2008]), which is the topic of this paper. 

In section two we describe how warm dark matter suppresses structure on small scales 
through free-streaming. We go on to discuss and compare two different methods of calculat- 

1 http://pan-st arrs.ifa.hawaii.edu/ 
2 https: / / www.darkenergysurvey.org/ 
3 Miyazaki et al. [2006] 
4 http://www. lsst.org/lsst/ 
5 http://sci.esa.int/euclid 
6 http://wfirst. gsfc.nasa.gov/ 
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ing the effects of non-linear evolution of the matter density field in the WDM case. We do 
not try to develop a new model for non-linear WDM structure as this is beyond the scope 
of this project. In the third section we calculate the weak lensing power spectrum for CDM 
and different models of WDM. In section four we report our predictions for a Euclid-like 
survey and the potential of such a survey to constrain the WDM particle mass. We take a 
flat ACDM or AWDM Universe throughout. Our parameters have the values of WMAP7 
(Larson et al. [2010]): the dimensionless Hubble expansion rate parameter, h = 0.71, the 
present day baryon density in units of the critical density, fib = 0.045, the present day dark 
matter density in units of the critical density, $7 m = 0.27, the cosmological constant energy 
density in units of the critical density, Qa = 0.73, the scalar primordial power spectrum 
slope, n s = 0.96, and the present day linear theory root-mean-square density fluctuation 
in 8 h^ 1 Mpc spheres, as = 0.80. 

2. The effect of WDM on the distribution of matter 

We begin by describing in more detail the evaporation of low mass halos through free- 
streaming and write down an approximate halo mass limit. We then go on to how we 
implement the halo model of the Universe in WDM scenarios. Finally we extrapolate the 
fitting formulae to CDM simulations in order to describe non-linear WDM structure. 

2.1 Free-streaming of WDM particles 

Free-streaming results in suppressed formation of dark matter overdensities smaller than 
the free-streaming length. We define a "free-streaming mass", Mf s as the mass contained 
within a radius of Af s with the average density equal to the background matter density. 
In other words, it is the mass of a dark matter halo that formed from an initial linear 
overdensity of size equal to the distance of warm dark matter free-streaming, Af s . We 
calculate the free-streaming length as (Zentner and Bullock [2003]): 

,21 1/3 r _ n _ 4/3 

Mpc (2.1) 



A fe ~ 0.11 



^wdm/i 



0.15 



WWDM 

L keV 



and define a corresponding 'free-streaming halo mass' as 

4tt /A fs x3 



where p m (z) is the background matter density at redshift z. 

The power spectrum of matter distribution is modified by this erasure of small mass 
halos. We modify the CDM linear theory power spectrum shape to the WDM case through 
multiplication by a transfer function acting on the linear theory density power spectrum 
(see Boyanovsky et al. [2008b] for a calculation of the transfer function for a general initial 
thermal distribution of DM particles - cold WIMP dark matter, thermal fermionic or 
bosonic dark matter): 

t <*Kct)' /2 -i i+ <^- 5/ " • < 2 - 3 > 



- 4 - 






250 eV 
500 eV 

1 keV . 
CDM 






\>» 

\ $ 

> v v 
1 \ 





10 10 10 10 10 

k, wavenumber [h Mpc~ 1 ] 



10 10 10 

k, wavenumber [h Mpc~ 1 ] 



Figure 1: Matter power spectra today (zero redshift), suppressed by free-streaming of WDM 
particles. The dashed lines are the linear power spectra corresponding to, from left to right in each 
panel, 250 eV, 500 eV, 1 keV warm dark matter and cold dark matter. The dashed lines are the 
same in both plots. The solid lines are the corresponding non-linear power spectra. The non-linear 
power spectrum in the left panel uses the halo model and that in the right panel was calculated 
using the haloht formula. The slight excess of power on large scales in the halo model comes from 
neglecting the effect of 'halo exclusion' (Smith et al. [2010]). 



(Bode et al. [2001]) with adjustable parameter \x = 1.12 (Viel et al. [2005]) found by fitting 
to the result of Boltzman integration (using CAMB by Lewis et al. [2000] or CMBFAST by 
Seljak and Zaldarriaga [1996]). The break-scale parameter, a in the case of thermal relic 
DM is 

■^(^"(wR^R" 1 "* <-» 

where S^wdm is the WDM density parameter, which in this paper equals the total dark 
matter density parameter, Odm> since all dark matter is assumed to be warm. The dashed 
lines in Fig. 1 show linear theory matter power spectra for a range of WDM masses. 
The lightest WDM particle mass shown (250 eV) causes the linear theory matter power 
spectrum to fall dramatically at a wavenumber significantly above 1 h Mpc -1 . As also 
shown in Smith et al. [2011], the matter power spectrum of WDM starts to turn off well 
above the free-streaming scale, which changes the slope of the power spectrum to fall much 
more steeply than n e g = log P(k)/ logk = —3, which is the slope for standard, bottom-up 
structure formation (White and Croft [2000], Knebe et al. [2003]). When, at higher k, 
the slope falls more steeply than that, bottom-up structure formation halts and top-down 
fragmentation might occur. This effect significantly alters the mass functions and since 
the effect starts at scales larger than the free-streaming scale, this influences halo masses 
above the 'free-streaming' mass (see Fig. 2). 

2.2 The halo model approach 

The halo model assumes the Universe is made up of halos with positions sampled from the 
linear theory matter distribution. As a result, there are two main contributions to the non- 
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Figure 2: The top figure shows Shcth-Tormen mass functions suppressed by free-streaming of 
WDM particles for two different warm dark matter scenarios: 1 keV particles (dashed red line), 5 
keV particles (solid green line) and CDM (solid black line) . The bottom figure shows the fractional 
difference between the WDM and CDM mass functions for the 2 different warm scenarios. The 
circles mark the halo mass on the x-axis that corresponds to two-times the free-streaming radius of 
each WDM particle (or eight-times the "free-streaming mass"). Data points from the simulation 
by Zavala et al. [2009] are marked with crosses on the upper panel: the red crosses correspond to 
the same WDM scenario as the red dashed line (only data points above the limiting mass of the 
simulation are shown). 

linear matter power spectrum (Ma and Fry [2000], Seljak [2000], Peacock and Smith [2000]). 
Firstly, the two-halo term, which dominates on large scales, encodes the correlation between 
different haloes and is equal to the linear matter power spectrum on large scales. Secondly, 
the one-halo term refers to the correlations within a halo and therefore depends on the 
density profile of the halo. Both terms depend on the number of halos as a function of halo 
mass, which can be found to a reasonable approximation using analytic arguments (Press 
& Schechter 1975) or more usually measured from numerical simulations and modified 
according to the underlying matter power spectrum (e.g. Bertschinger [1985], Sheth and 
Tormen [1999] and in the context of weak lensing, Cooray and Hu [2001] and Takada and 
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Jain [2003]). For reasons of simplicity we assume spherical collapse in our halo model 
calculations. 

We extend the halo model to WDM scenarios by modifying the ingredients. We use 
the WDM linear power spectrum described in the previous subsection to calculate a new 
mass function using the Sheth and Tormen [1999] prescription. We make the conservative 
assumption that the halo profiles are unchanged relative to CDM. It is in the one-halo term 
of the power spectrum that the effects of free-streaming of WDM are seen most strongly. 
This is because of the difference in the root-mean-square fluctuation, cr(R), which becomes 
suppressed at small R in a WDM Universe. 

For a detailed description of these calculations, see Appendix A. We choose not to 
normalise the mass functions because the result of free-streaming is the suppression of the 
formation of small structures, not by accreting the excess matter onto larger haloes, but 
rather by resulting in a larger uncollapsed background matter density. In contrast, we do 
normalise the halo bias by requiring 



This ensures that our model reproduces the linear predictions in the limit of k — > 0. 

We explore the effect of WDM on the mass functions in Fig. 2. As expected, the 
number density for the smallest haloes is reduced in the case of WDM. This is shown 
most visibly in the lower panel of Fig. 2 which shows the suppression on a linear scale. 
This is useful for comparison to the general assumption of the absence of haloes below the 
"free-streaming" mass. 

The points in the upper panel of Fig. 2 make a comparison to an N-body simulation by 
Zavala et al. [2009]. This indicates that the Sheth- Tormen mass functions underestimate 
the strength of this suppression. Indeed, comparing the simulated mass function for the 
1 keV particle to the Sheth- Tormen version, we find a very large discrepancy around the 
free-streaming mass. The simulated mass function declines to zero much more steeply 
than the Sheth- Tormen. This suggests that we underestimate the effect of WDM on the 
mass functions. New models need to be developed to calculate the WDM mass functions 
correctly, most likely on the basis of further WDM simulations. We have removed the first 
three points in the simulated mass function for the lowest three halo mass bins. This is 
because these are too affected by Poisson noise (Zavala et al. [2009]). 

2.3 Non-linear theory fitting formula (halofit) 

In the ACDM model, the non-linear evolution of the matter power spectrum can be cal- 
culated from the primordial linear matter power spectrum, using the halofit formula of 
Smith et al. [2003]. This method is a fitting function based on the halo model, but very 
accurately reproduces the results of many CDM numerical simulations. It is more accurate 
and more general than the fitting function from Peacock and Dodds [1996]. In this paper 
we make the large extrapolation of assuming that the halofit method is valid for WDM. In 
Fig. 1 we compare non-linear matter power spectra calculated using this method to power 
spectra calculated using the halo model formalism described in the previous section. 




(2.5) 
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Non-linear evolution of the matter power spectrum washes out the difference between 
WDM and CDM on small scales, which is apparent when comparing the linear spectra to 
the non-linear spectra in either panel of Fig. 1. This is due to the collapse of structures, 
which significantly boosts the power on small scales and more than compensates for any 
missing structure (Boehm et al. [2005]). The information contained in the matter power 
spectrum on cosmological parameters is encoded in the shape of the matter power spectrum. 
If the matter power spectrum is modified to smooth out the shape (e.g. by non-linear 
growth) then some of the information contained in the shape will be no longer be available, 
just as for any smoothing of an observable. In particular on small scales there is very little 
power to begin with, so when power cascades down from larger scales the little power is 
swamped. For any mildly noisy observation the information originally encoded in the small 
scales would be much harder to recover after the power had cascaded down, i.e. information 
would have been lost. 

3. Shear power spectra 

In this paper we consider how to measure the WDM particle mass using observations of 
cosmic shear power spectra. From an observer's point of view, the image of each galaxy in 
the Universe is distorted by gravitational lensing effects of all intervening matter. Therefore 
the cosmic shear power spectra are closely related to the matter power spectrum integrated 
over redshift. Future surveys are expected to use broadband photometry to estimate the 
redshifts of the observed galaxies. This should allow shear power spectra to be calculated 
at a range of different redshifts, and also allow cross power spectra between redshifts (see 
Csabai et al. [2003] for a review). We first calculate and compare cosmic shear power spectra 
calculated using the methods described in the previous section. We plot the comparison in 
Fig. 3. We then show power spectra for different WDM masses and compare visually with 
the uncertainties from an upcoming survey in Fig. 4. 

We consider a cosmic shear survey which has a number of galaxies per unit redshift 



(Smail et al. [1994]) where a = 2, /3 = 1.5 and z$ = z m /1.412, where z m = 0.9 is the median 
redshift of the survey. For example, this is reasonable for a Euclid-like survey (Amara and 
Refregier [2007]). We assume a photometric redshift uncertainty of 5(1 + z) = 0.05(1 + z) 
and take top-hat photometric redshift bins containing equal numbers of galaxies (as in 
Amara and Refregier [2007]). We use 35 galaxies per square arcminute and take / s k y , 
the fraction of the sky covered by the survey, to be a half. We assume there are no 
catastrophic outliers. We calculate power spectra at 20 spherical harmonic multipoles in 
the range 20 < t < 20000 (roughly 0.0065 < k < 6.5 at redshift 0.9, or down to scales of 
order 1 arcminute). 

The lensing power spectra are given by 



n(z) = z a e 
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Figure 3: In the upper panel we plot the weak lensing auto-correlation power spectra for tomo- 
graphic redshift bin 5, for 1 keV WDM particles and CDM. The lower panel shows the fractional 
differences between 1 keV WDM and CDM spectra for both methods: the thick solid (black) and 
the thick dashed (blue) lines were found using the WDM linear theory matter power spectrum with 
the non-linear contribution calculated with the halofit method and the halo model respectively. 
The thin dashed line was found using halofit with CDM in both panels. However it cannot be seen 
in the panel above, because it is only different from the thick (black) line by 1%. Note that the 
excess of power in the halo model on large scales is due to the normalisation of the halo bias and 
is unphysical. 



where xi{ z i) is the comoving distance to the lens at redshift z\ and W% is the lensing weight 
in the tomographic bin i: 



where 



Wi(zi) = p mfi 
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Figure 4: In the upper panel we plot the weak lensing auto-correlation power spectrum for tomo- 
graphic redshift bin 5. We show the power spectra for CDM (solid black) and three different WDM 
particle masses: 250 eV, 500 eV and 1 keV denoted with the dashed (red), dot-dashed (green) and 
solid (blue) lines respectively. The effect of the 1 keV particle is not visible on this plot and the blue 
line overlaps with the black. In the lower panel we plot fractional differences between the WDM 
and CDM power spectra shown in the upper panel. The (red) circles, the (green) squares and the 
(blue) diamonds represent the 250 eV, 500 eV and 1 keV particles respectively. Furthermore, here 
we plot the error bars for a Euclid-like survey in 20 multipole bins. 



and the subscripts s, 1 and Is denote the distance to the source, the distance to the lens 
and the distance between the lens and source respectively. 

In the halo model one can simplify the above equation using some approximations 
(e.g. Cooray et al. [2000]) in the halo model non- linear matter power spectrum so that the 
lensing power spectra are made up of one-halo and two-halo terms, 



c ij (i) = clM) + cfM) 



(3.5) 
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which are given by 
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Ti(l,z)= f """ M \h Z) MM, z)k.,(l, M, z)dM (3.8) 

(Cooray and Hu [2001]): where ki(l,M,z) is the 2- dimensional Fourier transform of the 
lensing convergence for which one needs to integrate the NFW halo density profile (Navarro 
et al. [1997]), truncating the integral at -Riso, the spherical radius, where the density 
contrast, A = 180. We use the analytic results given in Takada and Jain [2003]. The 
lensing convergence signal in redshift bin i is found by the following equation: 

K i (8,M,z l )= / ni (z s )K(0,M,z h z s )dz s . (3.9) 

We use the WDM linear theory matter power spectrum within the full halo model 
analysis, and calculate the resulting shear power spectrum, shown by the dashed line in 
Fig. 3 for the case of a 1 keV WDM particle mass, for the autocorrelation of the fifth redshift 
bin with itself. As expected, the cosmic shear power spectrum amplitude is reduced on 
medium to small scales due to the lack of contributions to the one-halo term from smaller 
halos. On mid-range scales the two-halo term is very slightly greater than the CDM 
power spectrum due to the renormalisation of the halo bias given in Eq. 2.5. The mass 
function is decreased on small scales by WDM and therefore the bias must increase on all 
scales to satisfy the bias normalisation constraint. For the matter power spectrum this 
exactly preserves the large scale power, but the lensing two-halo term weights the halos 
by their lensing effect (k in Eq. 3.8), which is largest for high mass halos, where the mass 
function was unchanged but the bias is larger. Consequently, as expected, the lensing 
power spectrum is increased on larger scales. 

We then use the linear theory WDM function within halofit to calculate the non- 
linear matter power spectrum, which we then insert into the standard cosmic shear power 
spectrum integrals. This is shown by the black solid line in Fig. 3. The effect is similar to 
that from the halo method but the deviation from CDM is about half the size. 

In using these two different methods we have attempted to somewhat span the range 
of possibilities given the present theoretical uncertainties. They all agree surprisingly well, 
to within a factor of two. To be conservative, we use the halofit method from now on, 
unless otherwise stated. 

For illustration we show the halofit results for three different WDM particle masses in 
Fig. 4. As expected, the suppression of the shear power spectrum is strongest and reaches 
the largest scales for the lowest DM particle masses. This reflects the scarcity of low mass 
lenses in the WDM Universe. The effect of WDM on the shear power spectrum (figure 
4) is more prominent than that on the non-linear matter power spectrum at redshift zero 
(figure 1). WDM has a greater effect on the linear theory matter power spectrum than 
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the non-linear matter power spectrum and therefore has a greater effect on the non-linear 
matter power spectrum at higher redshifts, where there is less non-linearity. Since the 
shear power spectrum is an integral of the matter power spectrum over redshift, WDM has 
a greater effect on the shear power spectrum than on the redshift zero non-linear matter 
power spectrum. We also see that the effect of halving the WDM particle mass is to more 
than double the perturbation to the CDM power spectrum i.e. the change in the deviation 
from CDM is non-linear in the inverse particle mass. We return to this in the following 
section. 

We compare the lensing power spectra in Fig. 4 to predicted statistical uncertainties on 
20 logarithmically spaced bins in I for this redshift bin. Here and in the following section, 
we assume a Euclid-like survey with 35 galaxies per square arcminute covering half the 
sky with shear measurement error on each galaxy of <r 7 = 0.35/\/2. (Refregier et al. 
[2010]). The comparison is only indicative, since there are strong correlations between I 
bins on small scales and between the 90 different redshift cross power-spectra (not shown). 
However, we already see that very high WDM particle masses around 250 eV will be easily 
ruled out. 

4. Forecasting constraints on the WDM particle mass 

In this section we forecast constraints on the WDM particle mass using a Fisher matrix 
(FM) analysis. We vary four other cosmological parameters in the cosmic shear analysis and 
compare constraints between the halofit and halo model approaches. We then combine the 
halofit results with forecast constraints from the Cosmic Microwave Background (CMB). 

We use CDM as our fiducial model and place limits on the inverse WDM particle mass 
so that our fidicual model has a parameter value of zero, corresponding to an effectively 
infinite WDM particle mass (in practice we use tbwdm = 100 MeV). We vary the following 
additional parameters: the total matter density today, O m , the primordial perturbation 
spectral index, n s , the primordial spectral amplitude, A s and the matter power spectrum 
shape parameter, V = e~ 2nhh Q m h. We use ^4 S instead of <r% to normalise the matter power 
spectrum for convenience when combining with the CMB constraints. For the CDM linear 
theory matter power spectrum we use a fitting formula for the matter power spectrum from 
Ma [1996] which uses the shape parameter T. Cosmic shear alone is very insensitive to f2b 
and h individually. 

The Fisher matrix for cosmic shear is given as: 



where Cij is the shear power spectrum of bin i correlated with bin j , (fi a are the parameters 
over which we marginalise and Cov(CV,(/), C mn (l)) is the covariance matrix of bins i, j, 
m and n. We use equation (14) in Takada and Jain [2004] to find the elements of the 
covariance matrix, after adding the shot noise term coming from the shear error, c 7 : 




(4.1) 



Cij (I) = C^ (I) + 5ij 



:2 

7 



(4.2) 
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Figure 5: 68 and 95 % confidence likelihood contours for both the halo model (blue dashed) and 
the halofit (green solid) approaches, marginalised over the power spectrum amplitude parameter 
A s and shape parameter T. The units on rn^h^M are e V -1 . 



where m is the average number density of galaxies in redshift bin i. 

The FM analysis assumes that the first non-vanishing order of the Tayor expansion 
of the log-likelihood is quadratic. However for CDM as the fiducial model, this is not the 
case for nt^ DM (or otwdm) as a parameter. Therefore the standard FM approach breaks 
down. As a result, the error found from FM analysis depends on the fiducial model chosen, 
errors are asymmetric around the fiducial value, and it turns out that the power spectra 
are completely insensitive to the inverse WDM particle mass at small values of m^y. 
Therefore it is not possible to use a very small step size in this parameter when calculating 
the derivative of the power spectrum in the Fisher matrix analysis. We discuss how we 
deal with this in Appendix B. Briefly, we tune the step size in the inverse WDM particle 
mass to equal the 1-sigma error bar on this parameter. 

We show the resulting parameter constraints on m\yDM) and n s in Fig. 5, marginalised 
over the power spectrum amplitude and shape. The halofit results are shown by the solid 
contours and the halo model by the dashed curves. The agreement between the uncertain- 
ties on each parameter is good, with the halo model giving slightly tighter constraints, as 
expected from the lensing power spectra shown in Fig. 3. The degeneracy direction between 
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Figure 6: 68 and 95 % confidence likelihood contours for the halofit approach using cosmic shear 
(thick green solid) and cosmic shear combined with Planck forecasts (red dashed line and thin black 
solid contours shaded red). The units on to^^ are eV _1 . 



m WDM an d is slightly different between the two methods. We have verified that on fix- 
ing all other cosmological parameters the degeneracy direction is extremely similar for the 
two methods, so the difference arises from the interplay between changing the WDM mass 
simultaneously with other parameters including the large scale amplitude parameter A s . 
The degeneracy direction between m^ DM and n s is as expected: increasing m^ DM causes 
a decrease in power on small scales, which can be partially compensated by increasing the 
power spectrum index n s . 

The marginalised errors that result from this analysis can be found in Table 1, using 
a one-tailed confidence limit for th-wdm- We ^ nc ^ approximately 30% smaller errors when 
using the halo model, compared to the halofit method. 

We include in the analysis the Fisher matrix for the Planck CMB experiment in ad- 
dition to the weak lensing Fisher matrix presented above. We derive the cosmological 
constraints from Planck following the description by the DETF Albrecht et al. [2006] (for 
a more detailed discussion see Rassat et al. [2008]). We assume that we will only use 
the 143 GHz channel data. This channel has a beam of #f w hm = 7.1' and sensitivities in 
temperature of T = 2.2 /iK/K and in polarization of P = 4.2 fiK/K. We take / s k y = 0.80 
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Fiducial value 


Halo model 
Shear only 


Ho 

Shear only 


lofit 

Shear + Planck 


^WDM 
™WDM [ eV ^] 


CDM 

io- 8 (« o) 


> 935 eV 
+0.00107 


> 645 eV 
±0.00155 


> 2500 eV 
±0.00040 


Q m 


0.27 


±0.0011 


±0.0016 


±0.0012 


n B 


0.96 


±0.031 


±0.035 


±0.002 


A s [h- 3 Mpc 3 ] 


1.3 x 10 5 


±0.119 x 10 5 


±0.183 x 10 5 


±0.005 x 10 5 


r 


0.18 


±0.0090 


±0.0123 


±0.0005 



Table 1: Table of 68 per cent confidence limits for a flat AWDM cosmology marginalised over 4 
other parameters in each case. The errors on the inverse particle mass, ??Wdmj are single tailed 
limits. 



as the sky fraction in order to account for galactic foregrounds. We use as a minimum I, 
^min = 30 in order to avoid problems with polarization foregrounds. As described in the 
DETF report (Albrecht et al. [2006]) we choose as fiducial parameter set: cu m , 9 S , lm4 s , 
Wb) ti b , t, where 6 S is the angular size of the sound horizon at last scattering, lnA s is the 
logarithm of the primordial amplitude of scalar perturbations and r is the optical depth 
due to reionization. After marginalization over the optical depth and the physical baryon 
density, ojb we then calculate the Planck CMB Fisher matrix in the parameters: T, A s , n s , 
by using the appropriate Jacobian of the involved parameter transformation (Rassat 
et al. [2008] or Eisenstein et al. [1998]). We point out that the Planck FM is computed 
without including the WDM paradigm. This is based on the assumption that the scales 
affected by WDM with m-wDM 100 eV are far away from the scales probed by CMB 
experiments. Therefore, while it adds quite strong constraints on the other cosmological 
parameters, especially on the curvature, the Planck FM does not add any constraints on 

m WDM' 

The parameter constraints are shown in Fig. 6 for the halofit method. We see that the 
very tight constraint on n s from the CMB breaks the degeneracy between n s and towdm 
from cosmic shear alone, to produce a tighter constraint on w-wdm- The single parameter 
marginalised results for halofit and the halo model methods are shown in Table 1. 

5. Discussion 

We have made the first estimate of constraints on the WDM particle mass from cosmic 
shear power spectra. We compared approaches to estimating the impact of warm dark 
matter on cosmic lensing. We found the halofit approach of Smith et al. [2003] more 
conservative that the halo model. We comment on the difference below. We find also that 
the effects on the mass functions (see Fig. 2) reach on mass scales significantly larger than 
the "free-streaming" mass. We made forecasts for parameter constraints from a Euclid-like 
cosmic shear survey combined with Planck and found that if the true cosmological model 
were CDM, we could place a limit of mwDM > 2.5 keV at 68 per cent confidence. 
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Our calculations are conservative in several ways. Firstly, we use a transfer function 
to relate the CDM linear theory matter power spectrum to the WDM linear theory power 
spectrum, from which we calculate the mass function. We showed a comparison with a 
mass function from simulations which suggest a much bigger supression of the number of 
lower mass halos. To improve our constraints it would be necessary to have a fitting formula 
to simulations giving the mass function as a function of WDM particle mass, extending 
to smaller masses than the present simulations. Ideally, a better physical model for the 
non-linear collapse of WDM structures would be developed. In addition, the new WDM 
models for the mass functions should be physical enough to translate also into the halofit 
method and explain the decrease in signal due to WDM in the non-linear power spectrum. 
Furthermore we assumed that the NFW profile is unaffected by the dark matter particle 
mass, whereas at the very least we expect the central core to be smoothed out. We assume 
that the effect of this on cosmic shear power spectra is subdominant to the decrease in the 
number of small halos. To overcome this limitation we would need to use a replacement 
for the NFW which gives the density profile as a function of WDM particle mass and halo 
mass. In order to truly account for the effects of WDM on halo substructure, one must 
run N-body simulations (see for e.g. Colin et al. [2000], Alam et al. [2002] and Zentner and 
Bullock [2003]). 

However, our calculations are optimistic in that we use mass functions and mass profiles 
obtained from simulations containing no baryons. To perform this measurement on real 
data it would be necessary to use results from simulations using baryons. These are in 
only the preliminary stages of development at the present time, and there is considerable 
uncertainty on the prescriptions used as well as on many free parameters that need to 
enter the calculations. The effect of baryons is to change the predicted cosmic shear power 
spectra, and to introduce uncertainty in the predictions. We assume for the purpose of 
this paper that these problems will be solved by the time we will be analysing a Euclid-like 
survey. We hope that studies like these will provide extra motivation to improve the quality 
of simulations using baryons. We also assume spherical collapse, which is less realistic in 
describing non-linear structure than ellipsoidal collapse, which has additional halo mass 
dependence in the collapse threshold. In other words, a further source of suppression of 
small scales is the decreased amount of collapse for small haloes, which tend to have higher 
ellipticity. 

Additionally, we normalise the halo bias in order to recover the linear matter power 
spectrum on very large scales. In fact, the correct approach would be to add a third term to 
the halo model, which would account for the uncollapsed matter. It would be necessary to 
specify the power spectrum and correlation of this uncollapsed matter, but it is not clear 
how best to determine these. We have chosen to normalise the mass- function weighted 
bias, following convention in CDM halo model calculations. However, the effect of WDM 
is to decrease the number of low mass halos and therefore the bias is increased on all scales 
to meet the normalisation condition. This causes an increase in power on mid-range scales. 
This may be the reason why the correlations between f2 m and rn^- DM are inconsistent be- 
tween the halo model and halofit (See Fig. 5). In order to correct this, further studies of 
WDM on halo bias should be conducted. These issues however only affect the halo model 
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calculations and do not influence our final constraint from the halofit method. A way to 
include the above suggested correction has been developed by Smith et al. [2011], where 
a term is added to the calculation to account for the contribution of uncollapsed matter, 
which solves the problem of bias normalisation. 

It is well known that Fisher matrices are an approximate method for forecasting un- 
certainties on parameters. The formal derivation of the Fisher matrix states that it gives a 
lower limit on the parameter errors. If the predicted data points are non-linearly related to 
the parameters, then we can expect likelihood contours in parameter space to be banana 
shaped, which is not captured within the Fisher matrix errors. Even when considering a 
single parameter the Fisher matrix approximation breaks down due to the flatness of the 
log-likelihood, giving a gradient of zero at the fiducial value of oi^dm = 0- We find this is 
a significant problem for the WDM particle mass. Small changes in WDM particle mass 
around our fidicual CDM model make negligible changes in the cosmic shear power spectra 
which makes it impossible to sensibly use a derivative taken with small step sizes in the 
Fisher matrix. We use CDM as our fiducial model in order to be able to compare our result 
to other methods (Lya, X-ray). We solve the above problem by using a larger step size, 
but the Fisher matrix method would ideally be replaced by a full Markov-Chain Monte 
Carlo exploration of parameter space. 

We have assumed that the CMB is insensitive to WDM when we calculated joint 
constraints with Planck, because the effects of WDM are only important on smaller scales. 
We have further assumed that only primary CMB anisotropies from Planck are used for 
constraining cosmology. However, gravitational lensing of the CMB itself can be used to 
reconstruct single-bin lensing tomography information on the lensing potential (e.g. see 
Lewis and Challinor [2006]). This would enhance the constraints on WDM further. 

Our cosmic shear calculations rest on a linear theory WDM power spectrum fitting 
formula which assumes dark matter particles are light and thermalised, like for example 
a light gravitino. This fitting formula can be straight-forwardly re-expressed for sterile 
neutrinos and our lower limit on the thermalised WDM particle can be converted to a 
lower bound on the sterile neutrino mass. For our final result using halofit plus Planck 
forecasts the limit is mwDM > 2.5 keV which can be translated to m us > 15.5 keV. If dark 
matter is made of sterile neutrinos created using the simplest mechanisms then constraints 
from cosmic shear can be compared to those from X-ray observations to test the consistency 
of the model. If cosmic shear were to place the upper limit we derive then it would be 
inconsistent with current X-ray lower limits and would rule out the simple sterile neutrinos 
as WDM candidates. 

Our forecast constraints from future cosmic shear studies are similar to the upper 
limits already being obtained from observations of the Lyman-a forest. The main advan- 
tage of using cosmic shear is that it does not rely on hydro-dynamical simulations of the 
intergalactic medium. However the Lyman-a observations are at higher redshift (around 
z ~ 3) than the cosmic shear (which roughly probes mass at z < 1). This is advantageous 
because the matter fluctuations are more in the linear regime and there is less reliance on 
simulations of non-linear theory. Also the linear theory matter power spectrum is much 
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more sensitive to the WDM particle mass. A deeper cosmic shear survey would gain from 
this effect. 

Cosmic shear power spectra are just one way to probe dark matter structure using 
gravitational lensing. However, higher order shear statistics such as the bispectrum are 
more sensitive to small scale clustering than cosmic shear and therefore would be expected 
to be more powerful in constraining the dark matter mass. Furthermore, higher order 
distortions of galaxies, such as flexion, also probe smaller scales and would be a better 
test. Bacon et al. [2009] find that galaxy-galaxy flexion is sensitive to structure of around 
10 9 solar masses. The methodology described in this paper could be extended to make 
predictions for these statistics. 
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A. Appendix: Non-linear power spectrum with the halo model 



We use the Sheth-Tormen (Sheth and Tormen [1999]) mass function (see also Takada and 
Jain [2003], Seljak [2000]): 



n(M)dM = [1 + {au)-P} /He"^ 2 ^ with v = ( 



D+(z)a(M) 

(A.l) 

where A = 0.353 as in Jenkins et al. [2001]. Also, a = 0.73 and p = 0.175. D + {z) above is 
the linear growth factor and a{M) is root-mean-square fluctuation in a sphere of volume 
containing a mass M in a smooth universe, calculated from the linear power spectrum. 
5 C is the critical threshold for a linear density fluctuation to collapse to form a halo (see 
equation A. 4), calculated from spherical collapse. 

From Takada and Jain [2003] the bias consistent with the Sheth-Tormen mass func- 
tions: 

av — 1 2p ( S c (z) \ 2 , , . 

b(M) = 1 + — — + — - -. . . with v = ^,. c ; L—r , A.2 

which must be normalised according to Seljak [2000]. 

The non-linear mass scale, M* (z) is defined to be the maximum mass that is undergoing 
non-linear collapse at a given redshift, which is given by the condition: 

^-{nTwk))^ 1 ■ (A ' 3) 

As mentioned above, we require the threshold for the linear density perturbations to 
collapse and form a halo. This was given by Henry [2000]: 



.1/3 

3(l2n) 2 / 3 l^mo -1 
5 c (n m , ,z)= V ; [1- 0.0123 log(l + x 3 )] with x=^ — -J- 

£Aj J. — i % 

(A.4) 

The density contrast of collapsed objects to the background mass density is taken to be 
A(flrn,o,z) = 180. Therefore the average halo density pt(z) = A(z)p m (z). 

The halo density profile is taken to be the Navarro-Frenk- White (NFW) profile with 
r s being the scale radius defined as the ratio of the virial radius and the halo concentration 
parameter: r s = ~, 

P(r) = -^—^ 7* ( A - 5 ) 

1 + 



c 



with 

M M,,) = « T . M . , (A.6) 

3 \ln(l + c(M.z))--£^ 



and halo mass is defined as 



M=^180p m (z) . (A.7) 
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The NFW profile, like other ingredients of the halo model is calculated from ACDM models. 
Other halo profiles have been explored by e.g. Alam et al. [2002]. The halo profile is a 
function of the halo concentration parameter, c(M,z) as in Cooray et al. [2000]: 



c(M,z) = a(z) 



M 



M*{z) 



b(z) 



(Ai 



where a(z) = 10.3(1 + z)~ - 3 and b(z) = 0.24(1 + z)-°- 3 . 

The resulting dimensionless power spectrum from Smith et al. [2003] (but also Cooray 
and Hu [2001], Seljak [2000] etc.), is given by: 



p#(M) 



l 



p2h 



(M 



(2tt) 3 

P\in(k,Z 



dM- 



dn 
d~M 



p(k,M, z) 

Pm,0 



1 2 



r dn , , , , N p(k, M, z) 
dM—b(M, z)^ — '-^ 
dM y 1 ' p mfi 



(A.9) 
(A.10) 



where p(k, M, z) is the 3D Fourier transform of the NFW halo density profile. Note that 
to get the full non-linear power spectrum we must add the 1-halo and 2-halo contributions: 
P ni (k, z) = P^(k, z) + P^(k, z). This is plotted in the right panel of Fig. 1. 



B. Appendix: Step sizes in the Fisher matrix for the inverse WDM par- 
ticle mass 

We first perform a simple likelihood analysis, only varying "T-\yDM> assuming all other 
parameters are fixed. We use this to illustrate our proposed method and assess how well 
it works in one dimension, where we can compare with the truth. 

As discussed in the main text, we choose to use the inverse WDM particle mass as the 
free parameter in our Fisher analysis because it has a finite value for our fiducial model, 
CDM. However as tends to zero the lensing power spectra become completely 

insensitive to the exact value, which makes it impossible to perform a traditional Fisher 
matrix analysis in which the step sizes used in calculating the derivatives are very small. 
This is illustrated in the left hand panel of Fig. 7 for the simple case where all other 
parameters are kept fixed. Clearly the choice of step size has a big effect on the result. 
We propose to use a step size such that it is equal to the error bar, as illustrated by the 
intersection of the dark (black) and light (cyan) lines in the left hand panel of the figure, 
shown by the (red) circle. Due to the convex nature of the function this cannot be done by 
iteration but instead is most simply done by making the figure we show. This can equally 
well be carried out in multi-dimensional space as in one dimension by working with the 
one-dimensional marginalised error bar. 

By working briefly in one dimension we can compare our result to the truth obtained 
by plotting the full likelihood as a function of m^py. This is shown in the right hand 
panel of Fig. 7. The solid line shows the exact likelihood and the dot-dashed line illustrates 
the Fisher matrix error obtained by setting the step size equal to the lcr error bar. The 
crosses and vertical lines show the 68% one-tailed upper limits from each curve. We see 
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Figure 7: Left: The la error calculated from Fisher matrices (solid black line), but with different 
differentiation step sizes used to calculate the gradient of the cosmic shear power spectra with 
respect to the parameter m^- DM . The red circle marks the point when the differentiation step size 
equals the la error on the inverse DM particle mass. We fix all other parameters. Right: The 
1-dimensional probability f (?n^yD]vr) ' s plotted with a solid black curve. The dot-dashed (blue) 
curve is the Gaussian probability distribution corresponding to the Fisher matrix error obtained 
when using a step size equal to the la error. The vertical dashed lines correspond to the 68.3 % of 
the probability. The red circle on the left corresponds to the vertical dashed blue line on the left. 



that the exact 68% upper limit is within a factor of two of the Fisher matrix limit using 
our proposal. Our proposal gives a more conservative result. We can also see from the 
difference in shape betwen the solid and dot-dashed lines quite how non-Gaussian the exact 
probability distribution is. We proceed by using this method in multi- dimensional space. 



